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Abstract 

Periodic orbit theory for classical hyperbolic system is very significant matter of how we can 
interpret spectral statistics in terms of semiclassical theory Although pruning is significant and 
generic property for almost all hyperbolic systems, pruning-proof property for the correlation 
among the periodic orbits which gains a resurgence of second term of the random matrix form 
factor remains open problem. In the light of the semiclassical form factor, our attention is paid to 
statistics for the pairs of periodic orbits. Also in the context of pruning, we investigated statistical 
properties of the "actual" periodic orbits in 4-disk billiard system. This analysis presents some 
universality for pair-orbits' statistics. That is, even if the pruning progresses, there remains the 
periodic peak structure in the statistics for periodic orbit pairs. From that property we claim that 
if the periodic peak structure contributes to the correlation, namely the off-diagonal part of the 
semiclassical form factor, then the correlation must remain while pruning progresse. 
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I. INTRODUCTION 



A. Chaotic property of classical systems and universality of spectral statistics 



Statistical properties of periodic orbits is one of fundamental problems in quantum chaos, 
i.e. quantum-classical correspondence for classically chaotic systems. Quantum chaos has 
been pursued in the last few decades □ QDDD . Especially in the context of semiclassical 



theory, the spectra 
correspondence 



statistics has been often referred to as an emblem of quantum-classical 
7|. According to the BGS ( Bohigas-Giannoni-Schmit ) conjecture 
there must be a universality of level statistics which goes along random matrix theory, that 
is, if the classical system has stochasticity then the spectral statistics for the corresponding 
quantum system obeys GOE one. But there is no self-evidence for illustrating the above 
facts. 



B. Periodic orbit theory and semiclassical analysis for spectral statistics 

Quantum - classical correspondence have been discussed mainly in spectral statistics for 
the last few decades. One of the objects for exploring the correspondence is density of states, 



d (E) = 6 ( E ~ E n) = (d (E)) + d osc (E) 



(1) 



where E is a energy and d (E), (d (E)) and d osc (E) are the density of states, the mean part 
and the fluctuating part, respectively. 



Before the terminology "Quantum Chaos" 



became used, Gutzwi 



trace 
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'ormula which was prologue to proceeding of periodic orbit analysis 
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ler derived the 
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dosc (E) = -LftJ^B^exp 

7T ft Kj 



-S, (E) 



(2) 



The trace formula represents fluctuating part d osc (E) of spectral density of states and it is 
derived by semiclassical evaluation of Feynman path integrals. 7, Sj, T 7 , B 11 k are index, 
action, period, amplitude factor and repetition number of periodic orbits, respectively. The 
amplitude factor B 1 consists of the Monodromy matrix M 7 and the Maslov index a 7 , 

1 



^\M,-I\ 
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7T 



exp ( i-<7 7 



(3) 



Admirably, the quantities which are factored in the trace formula are canonical invariant. 

One of the main objects of semiclassical analysis for spectral statistics is spectral form 
factor, 

K (r) = J°° ^|yy (d osc (E + I) d osc (E - I) ) e exp [-i2^{d {E)) W \ , (4) 

which was proposed by Berry and which is associated with semiclassical sum rule, i.e. 
bootstrap 

(d(E))=hmV2e(d 2 osc ^E)). (5) 



On the one hand, from the random matrices theory 



m 124, 

the spectral form factor is 



obtained by applying Fourier transformation to the two points level correlation. 

k goe = 2r - r In (1 + 2r) = 2r - 2r 2 + . . . , for r < 1, (6) 

provided r = T/T^, = 2nhd(E) is Heisenberg time. This is valid for the systems with 
time reversal symmetry. On the other hand, from semiclassical analysis, i.e. substituting 
the Gutzwiller's trace formula for the spectral form factor (0J), it becomes 



7,7 , K,K 



E,T 

(7) 

where (-)e,t means smoothing over energy E and time T. Employing Hannay - Ozorio de 
Almeida sum rule ( the principle of uniformity ) , Berry estimated diagonal approximation of 
the semiclassical form factor (JJJ). Then, for a chaotic system with time reversal symmetry, 
the diagonal approximation of the expression (JJJ) becomes 

K(t)~2t for t < 1. (8) 



Hannay - Ozorio de Almeida sum rule |21|, |22j represents bearing of periodic orbits' sum in 
hyperbolic system. 

TV — 1 - for T^oo. (9) 

^ I My - I\ 
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That is, in the context of the semiclassical theory flZJ), the counterpart of the first term of the 
random matrix form factor (jHJ) is interpreted as a result from the manner of the distribution 
of the periodic orbits. 
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C. Progress in analysis for correlation among classical orbits 



Contrary to the achievement above mentioned, there is the significant matter which had 
been left unanswered; in the semiclassical manner, as what we can take a counterpart of 
the reminder term of the random matrix form factor? As is mentioned in the preceding 
subsection, through the diagonal approximation of the semiclassical form factor, the first 
term of the random matrix form factor had been perceived as a result from the distribu- 
tional manner of the periodic orbits. However, the remainder terms of the random matrix 
form factor, — 2r 2 + • ■ ■ or — r log (1 + 2r), have not been barely accounted for by periodic 
orbits more than the fact that they can be described with the off-diagonal terms of the 
semiclassical form factor, i.e. combinations among classical periodic orbits except for its 
own combinations. Although it becomes clear that much more amount of off-diagonal terms 
than the diagonal terms reduce to the corrective effect through the cancelation in phases, 
we don't know how it works. In other words, the outstanding issue is how to explain the 
contribution which remain through the cancelation among the phases, especially in the con- 
text of classical dynamics. In other words, we have not known yet which combinations of 
periodic orbits make principal contribution to the off-diagonal part of the semiclassical form 
factor which corresponds with the remainder terms of the random matrix form factor. 

Toward understanding the relation between the off-diagonal terms of the semiclassical 
from factor (171) and the remainder terms of random matrix form factor © , many endeavours 

nnnn n 

have been made [2J, |24J, |2£|, I2i| • When it became the 21st century, Sieber and Richter |27| 
presented an approach: they assumed the pair of periodic orbits which one orbit had self- 
crossing and another enclosed one and the certain distribution which the periodic orbits 
obeyed, and using semiclassical analysis, they derived the counterpart of the reminder term 
of the random matrix form factor (JBJ). This pioneering approach and those following works 
share some ideas; the logarithm dependency of the target distribution, the uniqueness of 
partner orbits, etc.; and they prospered in derivation of second term of the form factor 



analytically by use of semiclassical theory [3, 2^, 0, l^ll 



Apart from the direct analysis for the spectral form factor, Shudo and Ikeda |36[ inves- 
tigated action pair difference in the kicked rotor system and they claimed significancy of 
the periodic peak structure of the action pair distribution in the context of dynamical lo- 
calisation. Concretely saying, they presented the classification of the pairs of (non-periodic) 
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orbits which contribute to the periodic peak structure with their symbol sequences: the 
orbits pairs have the symbol sequences which have period-4 cycle structure and have some 
other symbols inserted in the cycle structure. And they concluded that such classified orbits 
were the homoclinic orbits. 

Other approaches on correlation among periodic orbits are also interesting. For example, 



Efetof and Kogan 



derived the results of the random matrices theory via the non-linear 



;heir results is 
381 ] ). and the 



sigma model in disordered electronic systems. The diagonal approximation of 
equivalent to the lowest order perturbation theory ( Altshuler and Shklovskii 
higher order corrections would be equivalent to correlation which consists of pairs among 
periodic orbits. 



D. What is the object and what is the condition to be required for the periodic 
orbits' correlation? 

As just described, many works have been performed in this decade, but with all those 
works, it is not enough to understand the correlation. There remain several questions and 
rack of information to understand what is the correlation; it is not clear whether there is 
consistency between the periodic peak structure in the Shudo and Ikeda's work and the Sieber 
and Richter's idea; there is few researches to show whether or how "pruning" phenomena 
have an effect on the correlation; some researches posed the distribution of the periodic 
orbits which can contribute to the semiclassical form factor, but it is not clear that the other 
periodic orbits not to contribute to the semiclassical form factor can be vanished into nothing 
through cancelation among phases; there is a few works using actual numerical periodic 
orbits, but it is not sufficient to understand the correlation; if there is the correlation, then 
it becomes recognisable in statistics of not only implicit physical quantities, e.g. crossings 
of periodic orbits, but also obvious physical quantities, e.g. their actions and periods, but 
there is little direct analysis of the statistics for the obvious quantities. 

In accordance with the lessons above mentioned, we will discuss the object and the condi- 
tions for a correlation among periodic orbits as a rule of thumb. Requiring such conditions, 
we will find out the correlation in classically hyperbolic system. As already mentioned, sta- 
tistical property of periodic orbits is absolutely essential for understanding spectral statistics 
in semiclassical manner. However, it is too hard to solve the above problems simultaneously. 
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Hence, we must winepress into the following two points; (a) what is the variable for which 
the correlation must be appeared in the statistics?; (b) what is the condition which the 
correlation must satisfy? 



(a) what is the variable for which the correlation must be found in?: @ 

The statistics or the correlation to be observed must be measured in terms of the 
quantities which the semiclassical form factor contains explicitly, e.g. action S, period 
T, Monodromy matrix M, etc... . Otherwise, as a future task, it is hard to understand 
how the discovered correlation contribute to the semiclassical form factor and recreate 
the remainder terms of the random matrix form factor. From the viewpoint of explic- 
itly, for example, to treat crossing angle as the object for semiclassical analysis has 
disadvantage, because crossing angle is not contained explicitly in the semiclassical 
form factor and it is not canonically invariant and also it cannot be defined in some 
systems. And what is more, we believe that if the correlation which contributes to 
the semiclassical form factor does exist, then they must emerge also in an obvious 
expression, e.g. statistics among actions and periods for periodic orbits. Therefore, 
if we take a statistics about the semiclassical form factor among several variables, 
then we will give the highest priority to the action and the period. Because from 
the representation of the semiclassical form factor, it is clear that the most important 
quantities are not the stability and the Maslov index but the action and the period. 
Telling as it is, in the semiclassical form factor, there are the differences between the 
actions in the phases and the sums of the period in the delta functions. It is legitimate 
to consider that cancelation among phases and delta functions play a significant role 
for contributing to the semiclassical form factor. 

(b) what is the condition which the correlation must satisfy?: @ 

" property on the correlation what we must look for. In 



We impose "pruning-proo 
the context of "pruning" 39|, the correlation had not been thoroughly understood. 
Pruning is a phenomenon that structurally unstable periodic orbits are disappeared 
as a system parameter is var ying . And pruning is a general property for any hyper- 
bolic systems 

0,0,00 

LJ. W^. . The terminology "pruning" means, plainly 
speaking, breakdown of 1-to-l correspondence between periodic orbits and symbol se- 
quences, that is, number of periodic orbits decreases from how much ones "complete 



6 



system" should have essentially ( cf. subsections II B and C ). The word " pruning- 
proof property represents an immune property against such pruning process. It is 
an important matter whether such the correlation is independent of pruning. If a 
discovered correlation has not pruning-proof property, then such correlation cannot 
be called as universal property any longer, and then it is inadequate to adopt that 
correlation as significant contribution to the semiclassical form factor. Therefore, it is 
indispensable to examine whether a discovered correlation has pruning-proof property. 
If a correlation among periodic orbits constructs significant contribution to the semi- 
classical form factor, then the correlation must have a pruning-proof property. That 
is, pruning-proof property is a necessary condition for the correlation which gives some 
contribution to the semiclassical form factor. 

E. Aims of this paper 

The aims of this work are the following two points; (a) to find out a correlation among 
periodic orbits which must be measured in terms of the quantities of which the semiclassical 
form factor is apparently consists; (b) to confirm that the observed correlation is proof 
against pruning. If the above aims are achieved, then it turns out that we get a candidate 
for the correlation which contribute to the semiclassical form factor. 

In addition, we must mention that we ease the condition along which a correlation among 
periodic orbits must go. Meanwhile, we don't mind whether the correlation recurrence the 
form of the remainder terms of the random matrix form factor. Actually, with a numerical 
method, it is too difficult to recreate the equivalent of the remainder terms of the random 
matrix form factor with the correlation as well as only to find out the correlation itself. And 
then we intend to postpone understanding how a periodic orbits' correlation contributes to 
the semiclassical form factor in this paper. Therefore we would mind only what kind of 
correlation is there in classically hyperbolic system. 

F. Construction of this article 

In Section II, we will illustrate the objective system, specifically 4-disk billiard. In Section 
III, we will state the methods of analysis, especially the pair statistics for periodic orbits. In 
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(a) (b) (c) (d) 

FIG. 1: 4-disk billiards : (a) parameter setting, (b) the numbered 4 disks, (c) when r = d/2 = 0.5, 
(c) when r = d = 1.0 

Section III and IV, we will state the statistics for the single periodic orbits to support the pair 
statistics. In Section V, we will present the results for the above pairs' and singles' analysis. 
In Section VI, we will give summary and discuss the relation between the correlation and 
the pruning. 



II. SYSTEM 



A. 4-disk billiard 

Here we illustrate the object system, 4-disk billiard system j^j]. Billiard system is a 
model for idealised motion of a billiard ball ( point particle ) j^J which goes straight ahead 
in a given domain and collides elastically on a boundary. First, the definition of this system 
is given for the configuration of the domain Q in which billiard ball runs ( FigO]): 

QeR 2 \v, v = v uv 1 uv 2 uv 3 , (10) 

Vi = {(x, y) eR 2 \ (x- x ci ) 2 + (y- y a ) 2 < r 2 }, for z = 0,l,2,3, (11) 

where the radius r is common to all the four disks and 2d is the side of the square that each 
centres of disks is located on the four vertices respectively. Let us name each of the four 
disks as "i-th disk" in clockwise order from the top right disk ( Fig^(b) )• Each i-th disk is 
arranged its centre on the four vertices of the square respectively, that is, (x c0 , y c o) = (d, d), 
(xci,y c i) = (d, -d), (x c2 ,y c2 ) = {-d, -d), (x c3 ,y c3 ) = (-d,d). V { is the interior of each i-th 
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disk and T> is the total interior area of all the four disks. The domain Q is defined as the 
complement of R 2 with respect to T>. The billiard ball goes freely within the domain Q and 
collides elastically on the boundary dQ. When r = 1.0, the four disks touches to each other, 
and when r > 1.0, the four disks overlap each other. In the case of r > 1.0, the four disks 
divides the 2-dimensional configuration space into the two domains, one is enclosed by them 
and another encloses them. In any case, the periodic orbits can be in the area enclosed by 
4-disks, (-d <x,y,<d)U Q. 

In this system, it is easy to calculate periodic orbits in the sense of adequate correspon- 
dence between periodic orbits and symbol sequences 3|. For example, the billiard 
ball bounce with a succession of the i-th disks, in such way, "0", "1" , "2", "3", and then 
corresponding symbol sequence is "0123" ( see Fig|T]and Fig|2]). 



B. Pruning 

Pruning is, as its name suggests, trimming in symbol sequence's space. Primarily, in 
generally hyperbolic systems, any periodic orbit can be transliterate into an appropriate 
symbol sequence, but it is not necessarily vice versa. That is, it is not always true that 
there is periodic orbit corresponding to any symbol sequence. We describe this situation as 
"pruned" state or state that system has pruning property. In other words, pruned state is 
that there is physically impossible periodic orbit which is hypostatised from certain symbol 
sequence ( Fig|21). To put in another way, if there is a complete symbol sequence's space, 
we can regard it as a "tree" and call a symbol sequence as a "branch" of the tree. Pruning 
is defection of appropriate branches among the tree. 



C. Pruning in 4-disk billiard system 

Let us review the condition and the circumstance of pruning for the 4-disk billiards. All 
the four disks have the same radius r which is the system parameter to be varied from 
r = 0.5 to r = 1.0 simultaneously. The fact to be aware of is that variable is radius r, not 
distance 2d , apart from many of foregoing works. 

First, we consider not prunend state, i.e. corresponding symbol sequence's space is com- 
plete one. When r = 0.5, there is no eclipse, more precisely, looking from one disk toward 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

1 i 1 1 1 1 1 i 1 1 1 1 1 i 1 1 1 1 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

FIG. 2: To be pruned, i.e. to be bifurcated, orbits in 4-disk billiard: from left to right, 
orbit has '010212'/ '0101212' / '01012121' symbol sequence and has 6 / 7 / 8 collisions, from top 
to bottom, r = 0.9, 0.95, 1.0. They become pruned with going to downward in these figures, while 
they become bifurcated with going to upward. 



diagonal one, there is no obstacle to interrupt the view. Even on the morrow of the eclipse 
(r e = d/y2 = 0.70710678..), the corresponding space of symbol sequences has no pruning. 
Eventually up to r = r c = 0.90715514863324558..., pruning does not occur. 

Next, when r = r c = 0.90715514863324558..., first-pruned orbit touches on the 0-th disk 



42j ( FigEl ) • The periodic orbit to be pruned first among the whole set of periodic orbits 



has the symbol sequence, e.g. "0(30) n 31(01) n 0(10) n 13(03) n " and "3(03) n (10) n l(01) n (30) r 
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etc. (n — > oo), that is, after infinite number of collisions between 3rd-disk and 0-th disk, 
there are also infinite ones between the 0-th disk and the 1-st disk and again. 




0.5 



-0.5 




FIG. 3: Instant of first pruning in 4-disk billiard (r = r c = 0.90715514863324558...) 



Finally, after the fist pruning, the above first-pruned orbit becomes physically impossible. 
In this way, from r = r c , pruning occurs successively as r is increased ( FigJ21). 

In the 4-disk billiard system, pruning can be interpreted as reverse to bifurcations. That 
is, reducing the radius, physically impossible orbit come to realistic one, which we call this 
as "to be pruned" orbit, and simultaneously it becomes bifurcated. Conversely, as enlarging 
the radius, bifurcated orbits are integrated into a single orbit, and immediately such an orbit 
disappears as a ghost one, i.e. physically impossible orbit ( Fig 121). 



D. Numerical calculation of periodic orbits 

With the aid of not true periodic orbits but non-periodic orbits, which must be near 
to "true" one, the studies of semiclassical analysis of the spectral form factor have been 
achieved. We should not overlook that there might be "pruning" among the periodic orbits 
which has the same symbol sequence as that of the non-periodic orbit which is expected to 
contribute to form factor. Therefore, "actual" periodic orbits is indispensable to semiclassical 
analysis with the aid of the trace formula ((21). 

We calculated periodic orbits numerically, as the preparatory to the investigation of 



periodic orbits' statistics. Actual method to calcu 



ate any periodic orbit is achieved by min- 



imisation of whole ad referendum length of one |49|. |5l| . [52]. We perform this optimisation 
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by use of Newton-Raphson's method. To avoid 'chaotic behavior' of Newton-Raphson's 
method, we also use bi-sectional method before using the former method. First, making 
use of bi-sectional method, we can see that the numerical periodic orbit converges in the 
neighbourhood of the "true" periodic orbit within adequate range. Next, using Newton- 
Raphson's method, we converge numerical periodic orbits to the "true" ones within more 
narrow error range. In this way, we get "actual" periodic orbits numerically. 

Actually, the amount of the periodic orbits which we have calculated in 4-disk billiard 
systems are approximately 1,000,000 for every value r of radius which we taken. The limi- 
tation of this calculation is the frequency of collisions which billiard ball impacts against 4 
disks. We investigate through 20 time collisions. For generating long symbol sequence which 
produce appropriate periodic orbit, we treated long integer with the aid of CLN ( Class Li- 
brary for Numbers ) library j^]. We actually use periodic orbits for pairs' statistics to 17 
times collision for the reason why there are both limitations, the outperform of computers 
and the times. The number of the periodic orbits which has 17 times collisions is 949560 in 
r = 0.5 ( without pruning ), and 703999 in r = 1.0 ( with pruning ). 

III. NUMERICAL ANALYSIS I : STATISTICS OF PAIRS OF PERIODIC OBITS 
A. Analysis to do 

In this section, we will define statistics to analyse for correlation among pairs of periodic 
orbits. And we will state results of these numerical analysis in Section V. As we stated above, 
statistics of periodic orbits, specifically statistics of pairs among them, has became very 
significant matter in semiclassical analysis for spectral statistics. From consistency among 
Eqs. (JZj),©? and (|8ji . it is clear that a statistic on analysis beyond diagonal approximation 
are inevitable. Therefore we will investigate statistical properties of periodic orbits, more 
precisely, of pairs among them. In the semiclassical representation of the spectral form factor 
(J7J), there are the difference in action between pair of the periodic orbits ( 7, 7' ) included in 
the phase factor, S y — Sy , and the sum in period between pair of the periodic orbits included 
in the delta function, T 7 + Ty. Then, in working up statistics in reference to pairs among 
periodic orbits, we take particular note of not only differences between actions, S 7 — Sy, 
but also sums of periods, T 7 + Ty, in distinction from the former analysis. 
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B. The spectral form factor in billiard system 

As noted above, in the semiclassical spectral form factor (JJJ), differences between actions 
S 7 — Sy and sums of periods T 7 + Ty play significant role. Here we will illustrate that the 
matter of the correlation among periodic orbits, i.e. the contribution of the off-diagonal part 
toward the semiclassical form factor, reduce to that of the pair length statistics in a billiard 
system. For this end, we begin with the spectral form factor in a billiard system without 
external field. In a planer billiard system, utilising length of periodic orbit L 7 , action and 
period for any periodic orbit become the following: 

S 7 (E) = V2mEL 1 -> V2EL, p 
T 7 (E) = ^=-^=, 

where we substitute unity for mass of a billiard ball, m — 1. Then the semiclassical form 
factor is represented as 

4 ( L - ^) 

(12) 

that is, the spectral form factor becomes a function of periodic orbit's length L. Eventually 
question as to action difference and period sum are rehashed into that as to length difference 
and length sum. Therefore, statistics on length difference and length sum over periodic orbit 
is indispensable to consider a candidate for the correlation contributing to the semiclassical 
form factor. 

C. Distribution function of length difference / sum for periodic orbit pairs 

C AL±,m ( L ±) 

For the reasons which have been mentioned above, we define the statistics with regard 
to the pair statistics for periodic orbits. From brief consideration, the most fundamental 
statistics for pairs of the periodic orbits are distribution of length difference or sum for every 
two periodic orbits. 
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Here we define the following statistics: 

CAL„, m (L-) = Ul ± i I ™AL_ < L 7 - Ly < (m + 1) AL_}, (13) 
CAL + , m (£+) = ft{7 7' I mAL + < L 7 + L y < (m + 1) AL + }, (14) 

L± = L 7 ±L y , AL± = sup{L 7 ± Ly}/M, m = 0, 1, . . . , M - 1 e N. (15) 

7 

C^ L _ m (L_) / C^ L+ m (L + ) is the distribution function of length difference / sum for periodic 
orbit pairs. The distribution function of length difference / sum is the histogram which 
represents how many pairs of periodic orbits with their length difference / sum L±, are 
there in each divided segments, [mAL±, (m+ 1) AL±\. The domain of the histograms is 
from to supremum of L 7 ± Ly among given periodic orbits, i.e. [0, sup 7 {L 7 ± Ly}]. M is 
the number of bins which divides the domain of the distribution function, and m is a index 
to point out a bin. 7, 7' are the indices of periodic orbits which includes repeated ones. 
AL_ / AL + is a class interval which classify any periodic orbits pair according to their pair 
length L± and they are each derived by dividing the supremum of the length difference / 
sum among the periodic orbit pairs into the number of bins. In the remainder of this paper, 
we abbreviate C^ L _ (L_) / C*AL + ,m (-^+) as C~ {L-) / C + (L + ) respectively. 

Actually, we take the distribution C~ (L-) / C + (L + ) over the all periodic orbits with 
n — 17 collisions and we take the distribution for each the radius of the four disks r, from 
0.5 to 1.0 in the matter of 0.05. Each of the domains [0, sup 7 {L 7 ± Ly}] is segmentalised 
into M = 2 19 intervals and these distributions are not normalised. 

D. Joint distribution function of length difference and sum for periodic orbit pairs 

C AL ± ,m-,m+ ( L -' L +) 

In addition to the above statistics, there is need for more adequate analysis which consider 
length difference and its sum simultaneously. If the correlation which contributes to the off- 
diagonal part of the semiclassical form factor (J2J) is originated from the cancellation in the 
phase factor, Therefore, it is needed to take statistics for periodic orbit pairs which sums of 
their two periods are the same value. In the words of (|T2j) . we have to investigate statistics 
for pairs which sums of both lengths are the same value. 
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Therefore, we define the following joint distribution function : 

CL ± , m _, m+ (L-, L + ) = (t{ 7 ^ y | m_AI_ < L 7 - L y < (m_ + 1) AL_ , 

m + AL + <L 7 + Ly < (m + + 1) AL + }, (16) 

L± = L 7 ± Ly, AL ± = sup{L 7 ± Ly}/M±, m ± = 0, 1, . . . , M± — 1 e N. 

7 

We describe this function C^ L± m _ m+ (L_, L + ) as the joint distribution of length dif- 
ference and sum for periodic orbit pairs. In the remainder of the paper, we abbrevi- 
ate C^ i±im _ !m+ as C ± (L_,L + ). The joint distribution function, assemblage of 
both action differences and period sums, is two-dimensional histogram and represents how 
many pairs of periodic orbits are there in each divided cell, [m+AL + , (m+ + 1) AL+] x 
[m_AL_, (m_ + 1) AL_]. Each of the two ranges of the variables L± of the joint distribu- 
tion is from to supremum of L 7 ± Ly among periodic orbit pairs, [0, sup 7 {L 7 ± Ly}]. 

Let us explain the other factors of this joint distribution function. M_ / M + is the 
number of bins which divide length difference / sum for periodic orbit pairs respectively, 
and m_ / m + is the indices of each numbered intervals [m±AL±, {m± + 1) AL±]. 7,7' are 
the indices of periodic orbits which includes repeated ones. AL± is the class interval which 
categorise the pairs of periodic orbits with respect to, and it is each defined as the supremum 
of the length difference / sum sup 7 {L 7 ± Ly} among the periodic orbit pairs divided by the 
bin number, as is the case with C~ (L_) and C + (L + ) ( Eqs.(JT3J), (fT^j) ). We introduce to 
the coordinate system ( see Fig. |U). First, let us consider the original coordinate, L 7 -Ly, 
system. If the original axes are L 7 and Ly, and then the new axes are L_ = L 7 — Ly and 
L + = L 7 + Ly which are rotated by f ( FigJH ) and this coordinate system is the object to 
be considered. 

The joint distribution functions C ± (L_,L + ) are represented as a density plot. Any 
density plot for logarithm of joint distribution log C 1 * 1 (L_ , L + ) are expressed as a colour map. 
Actually, the whole domain of the joint distribution [0, sup 7 {L 7 — Ly}] x [0, sup 7 {L 7 +L 7 /}] is 
partitioned into M_ x M + = 2 12 x 2 12 cells and these joint distributions are not normalised. 
For each logarithm of the joint distribution \ogC ± (L_, L + ), the common colour map is 
applied to each of the different range which is determined according to the maximum and the 
minimum values among the function log C L + ). As in the case of preceding subsection, 
we take the joint distribution over all of the periodic orbits with n — 17 collisions and we 
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(a) (b) 

FIG. 4: Joint distribution function of length difference and sum for periodic orbit 
pairs: (a) coordinate system, (b) slicing joint distribution of length difference and sum along 
vertical direction = const, and along horizontal direction L_ = const. 

take the joint distribution for each the radius of the four disks r, from 0.5 to 1.0 in the 
matter of 0.05. 

E. Slice of joint distribution function of length difference and sum 

For searching the correlation, it is insufficient only to analyse the foregoing density plots 
because in logarithm of the distribution function log C^ 1 L + ), the fluctuating part, which 
may have detailed structure, become more subtle. Then we take the joint distribution in the 
form of anti-logarithm, C (L_,L + ). Moreover, the most straightforward idea to analyse 
the joint distribution are to behold their cross-section. 

For this reason, we also use following slice of the joint distribution function of length 
difference and sum for periodic orbit pairs. Briefly, the target statistics comes to the 1- 
dimensional histogram. The slice of the joint distribution C (L_, L + ) is altogether sec- 
tion of it for the 2-dimensional histogram. In equation either m_ or m + is fixed 
and we get the section of the 2-dimensional histogram. Therefore, the resultant his- 
togram represents how many pairs of periodic orbits which have the determinate value 
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L± = L 7 ±Ly G [m±AL±, (m± + 1) AL±] are there in each divided domains L T = L^Ly G 
[m T AL T , (m T + 1) AL T ]. The lines to be slit up the joint distribution is perpendicular to 
the axis L + = L 7 + Ly, in other words, parallel to the axis L_ = L 7 — Ly ( see FigJH (b) 
). Through this statistics, we will investigate the correlation among the orbits which are 
selected by the delta function. 

We use the same coordinate system as used in the preceding density plots. The distinction 
is the fact that the each distributions (L_, L + = consts) are taken along the direction for 
the difference direction L_ with respect to each the fixed coordinates L + = L 7 + Ly. As in 
the above density plot, the whole domain is segmentalised into M_ x M + = 2 12 x 2 12 cells. 
Moreover, for looking effect of taking histogram width, we segmentalise the whole domain 
also into M_ x M+ = 2 11 x 2 11 cells. Then we can see the structures in slices more detail 
than in the density plot itself. Also in this calculation, we take the histograms over the all 
periodic orbits with n = 17 collisions and we takn the slices of the above joint distributions 
for each the radius of the four disks r, from 0.5 to 1.0 in the matter of 0.05. 

IV. NUMERICAL ANALYSIS II : STATISTICS OF SINGLE PERIODIC ORBITS 
A. Length distribution for single periodic orbits G/\L,m (L) 

As an assistant to understanding the preceding pair statistics C ± L + ), we need also 
to take some statistics for the single periodic orbits. And we will ask the question whether 
behaviour of histogram can be varied according to manner of dividing its domain into sev- 
eral intervals. Because of the brevity and the translucent correspondence to the foregoing 
distributions, we define the following distribution: 

G AL , m (L) = tf{7 | mAL < L 7 < (m + 1) AL}, (17) 
AL = sup{L 7 }/M, m = 0, 1, . . . , M - 1 G N. 

7 

7 is the index for a periodic orbit which includes repeated ones. M is the number of items 
which divide length domain [0, sup 7 {L 7 }] , and m is the indices of each numbered intervals. 
AL is the widths of each segmented intervals. This length histogram for single periodic 
orbits indicates how many periodic orbits are there in each intervals [mAL, (m+ 1)AL]. 
By use of this statistics, it can be seen how the histogram behaves according to the width 
AL. As long as there are no confusion in this paper, we abbreviate Gal,™, (L) as G (L). 
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We calculated the histogram for the length among the single periodic orbits with 17- 
collisions. In this analysis, we adopted some dividing numbers M for the above purpose. 

B. Cumulative number of length spectrum for single periodic orbits N (L) 

In the above statistics, there is bothering problem, how long widths do we adopt in those 
histograms. Therefore, we need for the statistics which does not have to consider the width 
of the histograms. One of the most simple statistics which does not take a form of histogram 
is a cumulative number density of length spectrum. 

The cumulative number of length spectrum for the single periodic orbits is defined as 
how many periodic orbits which length is below L: 

oo 

JV(L) = J]e(L-L 7 ). (18) 

7=0 

It is well known that this function has its asymptotic form for large L (54J], 

N (L) ~ exp {h top L) /L = N asympt (L) , (19) 

where h top is topological entropy. 

Our aim is not to take the asymptotic form but to reduce the deviation from its. In the 
form of the cumulative function on its own, it is hard to recognise the fluctuation part as 
deviation from its asymptotic form because the number of N (L) is more tremendous than 
that of Nfi uc t (L) as increasing L. As a result, we consider the following function: 

N fluct (L) = N(L)- N asympt (L) . (20) 

That is, subtracting the asymptotic part N asympt (L) from the whole cumulative number 
density N (L), we can get the fluctuating part Nf[ uct (L). 

Actual procedure is the following: First, we took the cumulative number distribution 
N (L) of the length spectrum. The distribution expressed as exp (ht op L) / L in the asymp- 
totic form. Next, we performed the least-squares fitting with the function Nfn (L) = 
A exp (BL) I L, where A and B are fitting parameters, in a reasonable local domain and 
we employ this fitted function as the asymptotic part N asympt (L) (fTHjl of the cumulative 
function. Finally, we subtract the fitted function Nf it (L) from the whole distribution TV (L), 
and then we get the fluctuating part Nji uct (L) (|2U|) of the distribution. 

18 



V. RESULTS 



A. Distribution function of length difference / sum for periodic orbit pairs 

C AL±,m ( L ±) 

We took the distributions of length difference / sum for pairs C~ (L_) / C + (L + ) ( Eq. lfTSj) 
/ Eq. ()14|) ) among the periodic orbits with n = 17 collisions for each radius r ( FigEJand 
FigEl)- The domains [0, sup 7 {L 7 ± Ly}] are segmentalised into M = 2 19 intervals. And 
FigEl and FiglHl are detailed figures centring around the crest of the distributions FigEl 
and Figj7| respectively. From the numerical analysis, the distributions of length difference 
/ sum seem to depend on radius of disks at a glance ( FigEl and FigEl)- Beholding the 
length pairs' distributions from r = 0.5 to r = 1.0, we can find out variations with its radii 
in those distributions. When r = 0.5, it is easy to see the periodic peak structure, and 
furthermore, we recognise self-similar structures therein, While when r = 1.0, it is scarcely 
able to recognise such the structures. But it is not clear whether such periodic peak structure 
and even more self-similar one of those distributions truly become disappeared as increasing 
r. We need more detailed analysis to establish the facts of the matter. 

B. Joint distribution function of length difference and sum for periodic orbit pairs 

Further analysis is to take the joint distribution of length difference and sum for periodic 
orbit pairs C ± (L_,L + ) ( Eq.([TB|) ). As the above calculation, among the periodic orbits 
with n — 17 collisions, we took the joint distributions for each radis r. The each domains 
[0, sup 7 {L 7 — Lry/}\ x [0, sup 7 {L 7 + Ly}] are partitioned into M_ x M + = 2 12 x 2 12 cells. 
Beholding the joint distributions from r = 0.5 to r = 1.0 ( FigEl), we can find out variations 
with its radius in those distributions. When 0.5 < r < 0.6, rhombus pattern can be seen, 
while when 0.7 < r, it is hard to recognise such structure at a glance. But from the enlarged 
figures, the detail structure show up ( FigEU]) out of relation to the values of radius r. 
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FIG. 5: Distribution of length difference for periodic orbit pairs C AL _ m (L _): In the top, 
from left to right, r = 0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 6: Detailed distribution of length difference for periodic orbit pairs C AL _ m (L _) : 
in the top, from left to right, r = 0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 7: Distribution of length sum for periodic orbit pairs C^ L+ m : In the top, from 
left to right, r = 0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 8: Detailed distribution of length sum for periodic orbit pairs C~^ L+ m (L + ) : In the 
top, from left to right, r = 0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 9: (Colour Online). Density plots for joint distribution of length difference and sum 
for periodic orbit pairs C^ L± m m+ (L_,L + ): In the top, from left to right, r = 0.5, 0.6, 0.7; 
in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 10: (Colour Online). Detailed density plots for joint distribution of length difference 
and sum for periodic orbit pairs C^ L± m m+ (L_,L + ): In the top, from left to right, r = 
0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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C. Slice of joint distribution function of length difference and sum 

In addition to the above analysis, we took the slice of the joint distributions along the 
lines where the sum of pair orbit's lengths is constant, L + = L 7 + Ly = const.. We took 
these histograms with no-logarithm plot for searching detail structure. Needles to say, as in 
the preceding joint distributions, we referred to the periodic orbits with n = 17 collisions for 
each radius r. The each domains [0, sup 7 {L 7 — Ly }] x [0, sup 7 {L 7 +Ly }] are partitioned into 
M_ x M + = 2 12 x 2 12 , 2 11 x 2 11 cells. From the results of the slices of the joint distributions 
( FiglTTl and FiglT2l ). we can see the variations with its radii r. When 0.5 < r < 0.6, the 
periodic peak structure with broad periods can be seen, while when 0.7 < r, it is hard to 
recognise such structure at a glance. But in detailed structure, slices of the joint distribution 
of length difference and sum have the periodic peak structures in irrelevance to radius of 
disks r. 

D. Length distribution for single periodic orbits Gal,™ {L 7 ) 

For more understanding to the foregoing pairs' distributions, among the single periodic 
orbits with 17 collision, we took the frequency distribution for their lengths ( Eq.(|18|) ) for the 
case of r = 1.0. And we adopt the numbers M of all intervals, M = 700, 1400, to see effect 
by width of the histogram. Then we can see the periodic peak structures in these histograms 
( FigEI] ) an d there is significant deviation from the corresponding gaussian distribution. 
The periodic peak structure become to be split as varying the number of intervals. But 
rough structures of them remain, which is indication that there must be some structure in 
essential distribution independently of the width of the histogram. 

E. Cumulative number of length spectrum for single periodic orbits N (L) 

Because we would look at the problem of the manner of taking width of histogram from 
the outside, we analysed the cumulative number density iV (L) of the length spectrum (fTSJ) 
in the case of r = 1.0. In doing so, we referred to the periodic orbits which has 2 < n < 17 
collisions. After the least-squares fitting, we subtracted the fitted function Nf it (L) from the 
whole cumulative number N (L) and we get the fluctuating part Nfi uct (L) of the distribution 
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fig. 11: Slices of density plots for joint distributions C^ L± m _ m+ (L _, L+): the each 
domains are partitioned into M_ x M + = 2 11 x 2 11 cells; In the top, from left to right, r = 
0.5, 0.6, 0.7; in the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 12: Slices of density plots for joint distribution C^ L± m _ m+ (L_, L + ): the each domains 
are partitioned into M_ x M + = 2 12 x 2 12 cells; In the top, from left to right, r = 0.5, 0.6, 0.7; in 
the bottom, from left to right, r = 0.8, 0.9, 1.0 
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FIG. 13: Length histograms for single periodic orbits G/s.L,m (L)'- the domain of the his- 
tograms is divided into (a) 700, (b)1400. 
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FIG. 14: Cumulative number of length spectrum: (a) numerical result for cumulative number 
N(L), (b) fluctuating part Nfi ct (L) 



(HOJ). FigHH (b) shows the fluctuating part Nfi uct (L) of the cumulative number density of 
the length spectrum. There is exactly the periodic peak structure and moreover seemed to 
be self-similar structure. 
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VI. SUMMARY AND DISCUSSION 



A. Summary 

Let us summarise the main points that have been mentioned above. Clarification of peri- 
odic orbits' correlation is indispensable for semiclassical reasoning of spectral fluctuation. In 
the 4-disk billiard system, we performed the statistical analysis for the pairs among periodic 
orbits and for the single ones. From the viewpoint of brevity, we presented the statistical 
analysis on the differences and the sums which appear in the semiclassical spectral form 
factor: for the periodic orbit pairs, both the distribution function of length difference / 
sum C~ (L_) / C + (L + ) and the joint distribution function of length difference and sum 
C ± L + ) and for the single periodic orbits, both the length distribution G (L) and the 
cumulative number of length spectrum N (L). In those analysis, we found some universal- 
ity of deviation from relevant Gaussian distribution, concretely saying, the periodic peak 
structure. 

B. Distribution function of length difference / sum for periodic orbit pairs 

As the first step, we took the distribution of length difference / sum for periodic orbit 
pairs C~ (L_) / C + (L + ) ( Figs J5l6l / Figs J7|8l ). From the consistency with the analysis 
made by Shudo and Ikeda, it was anticipated that there must be periodic peak structure 
in those distributions. These histograms were calculated for each the radius of the disks r, 
from 0.5 to 1.0 in the matter of 0.05. In brief tendency, the histograms present gaussian 
distribution. Focusing attention on the detail structures, For small r(0.5<r<0.7), there 
seemed to be the periodic peak structure and they present self-similar ones, while for larger 
r ( 0.8 < r < 1.0 ), harder we recognize such structures. But we have no assurance that the 
periodic peaks structure are certainly disappeared in the large r. 
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C. Joint distribution function of length difference and sum for periodic orbit pairs 

Inspired from the configuration of the semiclassical form factor (J7J), we analysed the joint 
distributions of length difference and sum for the periodic orbit pairs (|16J1 . These joint 
distributions are taken in the form of the density plots ( Figs J9llOI ). As is the case with the 
above analysis, these density plots were taken for each the radius of the disks r, from 0.5 
to 1.0 in the matter of 0.05. Focusing attention on the detail structures, while for small r ( 
0.5 < r < 0.7 ), there seemed to be the periodic peak structure and they present self-similar 
ones, Although for large r ( 0.8 < r < 1.0 ) in the whole density plots, we can hardly 
recognise such as the above structures. But from the detail density plots, evidently, we 
become aware of the periodic peak structure. 

D. Slice of joint distribution function of length difference and sum 

For intelligibleness, we taken the slices of the joint distributions of length difference 
and sum, i.e. the sections of the above density plots ( Figs ll 11121 ). And we taken these 
distributions in anti-logarithm for vertical axis of the histogram. Also in these sections, 
we can recognise that the periodic peak structures remain as varying r. Therefore, we have 
further confirmation that there is some universality among the above density plots. However, 
the manner of taking histogram's width is left open problem. 

E. Length distribution for single periodic orbits Gal,™ (Lj) 

For the understanding for the periodic peak structures as observed in the above analysis, 
we took the frequency distribution G (L) for the length of the orbits which has n = 17 strike 
points ( FiglTHl ). This analysis performed in the value of r = 1.0. Indeed, the periodic 
peak structure can be seen in this analysis. But the spread of the peak seemed to be varied 
according to how we take the widths of the histograms. If the periodic peak structure has 
self-similar structure, it is critical matter how to take the histogram width. 
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F. Cumulative number of length spectrum for single periodic orbits N (L) 

To avoid the bothering 'histogram width' problem in the anterior analysis, there is need 
for the analysis which is free from the histogram. We took the cumulative number for the 
length spectrum in r = 1.0. From it, we ultimately reduced the fluctuating part Nfi uct (L) 
through the least-square fitting. In the fluctuating part, we observed the peak structure 
and moreover recognised self-similar structure. Therefore, we assure there are indeed the 
periodic peak structure of the length histogram and besides that of the joint distributions. 



G. Picture of periodic peak structure 

From the above-mentioned analysis, we can claim as the following: In the joint distribu- 
tion of length difference and sum C ± L + ), there are in fact the periodic peak structure 
in a large, while in the histograms of length difference / sum C~ (L-) / C + we cannot 
say that the periodic peak structure exists irrespective of r. 

As the above mentioned, the periodic peak structure of the joint distribution has pruning- 
proof property. If such structure contribute to the off-diagonal part of the semiclassical form 
factor, the off-diagonal part has also pruning-proof property. This means the following: 
The correlation is the constraint condition for which preserve the periodic peak structure, 
However we do not know yet how many pairs of periodic orbits contribute to the periodic 
peak structure, and whether the number of periodic orbits which contribute to the form 
factor can be varied as success of pruning. Stated another way, we do not know more of the 
fact that the pruning fronts is a constraint condition for the correlation. 

We consider the fact that there is seemingly contradictory behaviour between the distri- 
bution of length difference / sum and the joint distribution of length difference and sum. 
From the density plots ( FigEU), it appears that the peak lines run along a oblique direction 
against the lines, L 1 — L 7 > = const, and L 7 +Ly = const.. Furthermore, there is an apparent 
relation between the distribution C~ (LJ) / C + (L + ) and the joint distribution L + ) 
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as follows: 

M+-1 

Cl L _, m _ {L-) = J2 C L ± , m ., m+ (L-, L + = m + AL + ) , (21) 

m+=0 
M_-l 

Cl L+ , m+ (L + ) = J2 C L ± ,m-,m + (L- = m_AL_, L + ) . (22) 

m_=0 

Consequently, the periodic peak structure in the joint distribution for the periodic orbit 
pairs seem to be cancelled when they are summed up to the distribution of length difference 
/ sum. From the perspective of the semclassical form factor (j?J), it is not so important to 
discuss absence or presence of the peak structure in the distribution of length difference / 
sum. 

In the second place, we consider the cause of the rhombus pattern in the density plots 
for the joint distributions of length difference and sum. And so let us review the coordinate 
system again ( see Fig0] (a) ). The rhombus pattern are formed along the lines L 7 — Ly and 
Lj + Ly. This fact indicates that the correlation ( periodic peak structure ) exists originally 
in the length histogram for single periodic orbits. 

From the last analysis, we got the periodic peak structure from the statistics which 
does not assume the form of a histogram. Therefore, we can claim that the preceding 
length distribution for single periodic orbits which is shown as a histograms have indeed 
the periodic peak structure. This indicates that such deviation from asymptotic behaviour 
exp (h top L) / L must contribute to the off-diagonal term of the semiclassical form factor. The 
statistics which we have analysed has probably self-similar structure, then we must treat 
the class interval of the histogram in a prudent manner. 

Toward further understanding, we must make the correlation's picture more detail. It is 
natural to question whether the periodic orbit pairs which contribute to the periodic peak 
structure can be changed as varying r. Owing to structural stability of the periodic peak 
structure, all we can say is that the combination of the periodic orbits which constructs par- 
ticular peak among periodic peak structure must be fixed. For this reason, such contributing 
orbits must be structurally stable. Against the interpretation mentioned at preceding sen- 
tences, structural unstable pairs, which consist of bifurcated orbits ( e.g. FigOl and FigEl 
), after all, do not contribute to particular peak. Because compared to length difference of 
structural stable ones, those of structurally unstable ones are variable and ultimately come 
to vanish as increasing r. 
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Next, we would like to mention a manner of sorting out the contributions to the periodic 
peak structure: which kind of set must we classify the contribution to such structure accord- 
ing to, a set of single periodic orbits or a set of pairs of periodic orbit ? From the histogram 
of length for single periodic orbits and the preceding considerations for it, we speculate that 
assemblage to be performed is to censor in light of a set of single periodic orbits. 

One of special properties of billiard system is the inseparably relation between action 
and period (fT2j) . Therefore, at a glance, it seems to be reasonable that the above statistical 
properties are highly special. However, essentially for a general system, action and period 
are closely related to each other. 

T(E) = ^§1. (23) 

Consequently, the preceding idea for the correlation, i.e. the periodic peak structure, must 
be extended. 



H. Conclusion 



We have performed numerical calculation of not non-periodic orbits but actual periodic 
orbits in the 4-disk billiard system. And we analysed the distributions of length difference 
and I or sum among the periodic orbit pairs and the histogram and the cumulative function 
of length spectrum for the single periodic orbits. From these analysis, it is clear that the 
periodic peak structures are observed in the joint distribution of length difference and sum, 
that is, the correlation to be required is observed in the statistics for the quantities which 
the semiclassical form factor explicitly consists of. And this is ultimately a result from 
the correlation which the length spectrum for single periodic orbits has. And notably such 
property remains unaffected by pruning. Furthermore, we ask whether pruning or bifurcation 
can become constraint condition for the correlation, and we conclude if either of the pair 
of orbits come into bifurcation then such the pair cannot contribute to the periodic peak 
structure. That is, if the periodic peak structure has very significant contribution to the 
form factor, we do not have to care bifurcation. 
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